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ABSTRACT 

We show that the distance measures (such as the luminosity and angular diameter distances) 
are linear functional of the equation of state function w(z) of the dark energy to a fair degree 
of accuracy in the regimes of interest. That is, the distance measures can be expressed as a 
sum of (i) a constant and (ii) an integral of a weighting function multiplied by the equation 
of state parameter w(z). The existence of such an accurate linear response approximation 
has several important implications: (a) Fitting a constant-ui model to the data drawn from an 
evolving model has a simple interpretation as a weighted average of w(z). (b) Any polynomial 
(or other expansion coefficients can also be expressed as weighted sums of the true w. (c) A 
replacement for the commonly used heuristic equation for the effective w, as determined by 
the CMB, can be derived and the result is found to be quite close to the heuristic expression 
commonly used, (d) The reconstruction of w(z) by Huterer et al. (2002) can be expressed 
as a matrix inversion. In each case the limitations of the linear response approximation are 
explored and found to be surprisingly small. 

Key words: cosmology:theory - methods: statistical -cosmological parameters. 



1 INTRODUCTION 



Current cosmological data suggests that the expansion of the uni- 
verse is accelerating (e.g. Perlmutter et al. 1999; Riess et al. 1998; 
Efstathiou et al. 2002; Lewis & Bridle 2002; Melchiorri et al. 
2002). A number of models have been proposed to explain this fact, 
the simplest of which is the cosmological constant. While a cosmo- 
logical constant is enough to explain the current data, its constancy 
leads to a fine tuning problem (Sahni & Starobinski 2000; Peebles 
& Ratra, 2002; Padmanabhan, 2002b). A simple phenomenological 
generalization of the cosmological constant is to model the dark en- 
ergy component that drives the acceleration as an ideal fluid with 
an equation of state given by P = wp in which the equation of 
state parameter w is allowed to vary with time. In this parameter- 
ization the cosmological constant corresponds to w = —1, while 
for other viable models —0.6 < w < — 1 at the present epoch. 
Such a parameterization indeed arises naturally in several mod- 
els, such as quintessence (Ratra & Peebles 1988; Wetterich 1988), 
K-essence (Armendariz-Picon et al. 1999) and a tachyonic scalar 
field (Gibbons 2002; Padmanabhan 2002a; Padmanabhan & Roy 
Choudhury 2002; Bagla et al. 2003). However the precise form of 
w{z) is model dependent. 

In view of this there has been considerable interest in attempts 
to summarize the current (and future) data in terms of a few num- 
bers. There are several ways of doing this, such as using a polyno- 
mial approximation for w(z) (e.g. Weller & Albrecht, 2002) or in 
terms of derivatives of the expansion factor (called Statefinders by 
Sahni et al. 2002). To be able to fit data using a polynomial form we 
have to choose a low order polynomial since realistic data contains 



only a finite amount of information. On the other hand we need 
to use enough polynomial orders so that the data is adequately de- 
scribed. Saini et al. (2003) investigate in some detail the issue of 
how to go about this and show that the current data and near fu- 
ture data seem to require at most a low order polynomial. (This, of 
course, does not imply that the true equation of state is also of a 
low order polynomial form.) 

These methods attempt to capture the effect of an unknown 
function w(z), which is equivalent to infinite number of parame- 
ters, in terms of a finite (and often small) number of parameters. In 
the limited redshift range probed by supernova data the true varia- 
tions in w(z) might be such that a low order polynomial fit is a rea- 
sonable approximation. Alternatively one can expand the function 
w{z) in terms of a some set of basis functions which are complete 
in the given redshift interval. If the basis functions are chosen judi- 
ciously only a small set of expansion coefficients will be required 
to describe the function w(z) within the limits of the experimen- 
tal accuracy. Owing to the random noise in cosmological data any 
of these choices could fit the data adequately in a maximum likeli- 
hood sense. Since the distance measures are integrals over nonlin- 
ear expressions involving w(z), any parameter determined by such 
a maximum likelihood analysis represents some non-trivial aver- 
age of the true w(z). Our main aim in this paper is to show that the 
actual functional relation of the cosmological distance measures to 
w(z) is not far from being linear. This allows us to deduce an ap- 
proximate relation between the fitted coefficients and the true w ( z ) . 
This allows one to obtain a quantitative description of the averaging 
involved in reducing a function w(z) to a finite set of numbers. 

The plan of this paper is as follows. In Section 2 we derive the 
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Figure 1. A surface plot of the kernel K w for supernovae distributed uni- 
formly up to z = 2.0. 

linear response approximation relating the coordinate distance and 
the equation of state of the dark energy, and show that within a rea- 
sonable range of parameters for the dark energy it works to better 
than ~ 2%. In Section 3 we use this approximation to relate the 
fitted coefficients of a polynomial form to the true equation of state 
for the case of fitting to the luminosity distance, and show that they 
are related through a weighted integral of the true equation of state. 
We then generalize this result to the case of an arbitrary functional 
form which is linear in the parameters. In Section 4 we derive a 
similar relationship for CMB where, in the simplest case, only one 
measurement of the angular distance to the last scattering surface is 
available. In Section 5 we show that the linear approximation also 
enables a non-parametric estimation of w(z). Our conclusions are 
presented in Section 6. 



2 RESPONSE OF GEOMETRY TO THE EQUATION OF 
STATE 

The luminosity and angular diameter distances depend on the dark 
energy through the coordinate distance r[w, z] according to the 
equations: Dl(z) = (l + z)r[w, z] and Da(z) — r[w, z]/(l + z). 
The square brackets in this notation explicitly shows that at any 
redshift z the coordinate distance requires full knowledge of the 
function w(z) In a flat universe the coordinate distance r[w, z] is 
given by 

r[w,z) = (1 + 3) 1/2 J Z dxx- 3/2 [g + Q[w,x]]- 1/2 (1) 
where x = 1 + z, g = fi m / (1 — Jim), 

rl+z 

Q[w, z] — exp[3 / dxw(x)/x], (2) 

and we have set c and Hq equal to unity. To explore the behaviour 
of r[w, z] when different equation of state functions w(z) are used, 
we need to understand the sensitivity of r to w(z). This can be 
characterized by the functional derivative of r with respect to w(z). 
Since this is not a routine weapon in the arsenal of the astronomer, 
we shall briefly introduce the concept before applying it. 

In the case of a real function f(x), the sensitivity of the func- 
tion to the independent variable x can be characterized by the 
derivative df/dx = f'(x). Broadly speaking, a large value for 
f'(a) indicates that / is relatively more sensitive to the indepen- 
dent variable around the point x = a; and a small value for /'(a) 



will indicate relative insensitivity of / to the independent variable 
around x = a. This follows directly from the definition of deriva- 
tive of a function 

f(a) = lim f(* + ^-f^ . (3 ) 

In the case of / depending not on a single variable but on a func- 
tion p(x) we need to study how / changes if the function p(x) is 
changed slightly around a point x = b. This is best done by chang- 
ing the function p(x) to a new function pi (x) = p(x) + e5i>(x — b) 
which adds a "spike" at x — b with a strength proportional to e. We 
can now evaluate the value of / for both p(x) and pi(x). The dif- 
ference in the numerical values of / in the limit of e — » is a 
good measure of the sensitivity of / to the functional form of p(x) 
around x — b. More formally, this functional derivative is defined 
by 

Sf[p,x] _ f[p + eSu{x -b),x\- f\p,x] 
Sp{b) ™ e ' { ' 

Just as the ordinary derivative of a function depends on the location 
at which it is evaluated, the functional derivative depends both on 
the form of p(x) around which the sensitivity is measured, x = a, 
as well as the point x = b at which the input function is perturbed. 

For the purposes of this paper we need the response of the co- 
ordinate distance r[w, z] at a redshift ztoa change in the equation 
of state at a different redshift, z'. This can be computed around a 
given fiducial w(z) = w Rd (z) from the functional derivative 

5r[w fid ,z] _ Um rlw^ + eSojz-z')^] -r[w Rd ,z] 
Sw(z') e^O e 

defined exactly as in equation J4j (also see Huterer & Turner, 2001 
for a similar application of this idea) 

For the rest of this section we switch from discussing the co- 
ordinate distance to the luminosity distance and change the inde- 
pendent variable from redshift z to x = 1 + z. Multiplying by 
(I + 2) we obtain the response function for the luminosity distance 
8Dl [w fid , x]/8w (x') = K w (x,x'). The subscript w on the ker- 
nel denotes that it depends on w Rd , the fiducial equation of state 
around which the approximation holds. Evaluating the expression 
in Eq|5]we obtain 

!_3£(i+£i^2 r x d v Ql™ lid ,y] 
2 *' Jx ' v 3/2 ( a +Q[^,y}) 3/2 . (6) 
(for x < x') 

The kernel K w is a function of two parameters, the redshift at 
which w(z) is perturbed and the redshift at which we consider the 
change in the luminosity distance. A surface plot of this function 
is shown in FigQ Since the effect of varying w(z) at a redshift 
z' is felt only at z > z',the kernel is zero in half the plane. For 
small Sw we can use this result to approximate the calculation of 
the luminosity distance as 

D L [w &d + 5w,z] ~ D L [w Rd ,z] + 6D L (7) 

= Dl [w fld , z] + / K w (z, z')8w{z')dz . 
Jo 

To be useful this linear response approximation should hold 
to a good accuracy for a reasonable range of w(z). The SuperNova 
Acceleration Probe (SNAP) survey is expected to observe about 
2000 Type la supernovae (SNe), up to a redshift z ~ 1.7, each 
year (Aldering et al. 2002). A single supernova will measure the 
luminosity distance with a relative error of ~ 7%. If we bin the su- 
pernovae in redshifts interval of 0.02, this will give a relative error 
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Figure 2. This Figure illustrates the accuracy of the linear response approx- 
imation by plotting the relative difference in the luminosity distance com- 
puted using the exact expression compared to that from the approximate ex- 
pression in Eq 7. This is calculated for the special class of models in which 
the equation of state is a linear function of redshift, w(z) = wq + w\z. The 
shaded regions give the relative error Sdi/di < 0.005,0.01,0.02,0.05 
from the inner hatched region diagonally outwards. The two straight lines 
mark the regions in which the — 1 < w(z) < in the range < z < 1.5. 



in the luminosity distance of about ~ 1%. Saini et al. (2002) show 
that given this level of precision and given the present uncertainty 
in the value of fl m , the data seems to require at the most a linear 
polynomial order in w(z). To show how accurately Eq[7]gives the 
luminosity distance in this restricted case, in Fig we plot in the 
wo-wi plane the percentage difference between the true luminos- 
ity distance and the one computed through Eq0at z — 1.5. The 
kernel K was computed for w (z) — —0.5 in this calculation. 
Within the region —1 < w(z) < the approximation is better 
than 2%. Therefore, the accuracy of the linear approximation is 
quite acceptable with respect to the projected uncertainties in the 
luminosity distance obtained from a SNAP class experiment. The 
approximation works even better at smaller redshifts but less so 
at higher redshifts, but the measurement errors are expected to be 
larger here. Therefore, the conclusions drawn from this approxima- 
tion are expected to remain valid until an unprecedented accuracy 
is achieved in the measurement of luminosity distance. 



3 INTERPRETATION OF THE FITTING PARAMETERS 

Since we have no fundamental understanding of the nature of dark 
energy, it is necessary to assume some suitable, versatile, functional 
form for w(z) to fit the cosmological data. Since the luminosity 
distance depends on the equation of state through an integral re- 
lation, and the cosmological data is often fitted through maximum 
likelihood method, we can always add a small fast varying term 
to any well fitting w(x), while still retaining a good fit , Due to 
this integral dependence any flexible parameterization eventually 
recovers only some integrated property of the underlying equation 
of state. The simple polynomial forms for w(z) usually assumed 



For example w(x) and w(x) + 9 s'm(x/9), where 9 0, would both 
fit the data equally well though none of the derivatives agree. Although this 
is a contrived example this point has been well illustrated more realistically 
in Maor et al. (2000). 



for the purposes of fitting the cosmological data are often viewed 
as the first few terms in a series expansion of w(z). Since the be- 
haviour of w(z) is not necessarily polynomial like, the recovered 
coefficients of the polynomial are not related in a simple manner 
to the true coefficients of the series expansion. In this Section we 
use the linear response approximation described above to compute 
the expectation value for the coefficients of the polynomial approx- 
imating the w(z), and relate them to the "true" input w(z). This 
will show that such fitting functions serve a useful purpose of mea- 
suring some integrated properties of the true equation of state of 
the dark energy. 



3.1 Interpreting the fit with a constant w 

The simplest fit to the cosmological data is the constant w model, 
w Rt (z) = Wq*. We are interested in relating this to the true 
w(z) — w tIUC (z). For an arbitrary equation of state this relation is 
non-trivial and non-linear but given the approximations considered 
in § |2| we can construct the relationship analytically. Let us first 
consider fitting the luminosity distance to redshift. This requires a 
knowledge of the present day matter density fi m as well and, for 
simplicity, we shall assume that this is known. (When Q m is un- 
known, both Q m and w become biased; see Maor et al. 2001). Our 
approach could be extended to the case in which f2 m is not known, 
however this is beyond the scope of the present paper. 

In the maximum likelihood reconstruction the quantity 



X 



dz 



{DT c (z)+n{z)) 
a(z) 



(8) 



is minimized, where n(z) is the noise in the measurement at red- 
shift 2 and u(z) is the variance. (We have replaced the conventional 
summation with an integral.) We now approximate the fitting func- 
tion and the luminosity distance by 



Df(z) ~D L [w &d ,z] + Aw 



D L {z) ~D L [w Rd ,z] + 



K w {z,z')dz (9) 
' K w {z,z")Sw(z")dz" (10) 



where we have assumed a constant w fld and defined 

w Rd and 5w{z) = w tIUC (z) — w fld . Minimizing \ 2 and taking the 

expectation value we obtain 



Aw 



where 



$ w (z')Sw(z')dz 



JJ K w (x, x')K w (x, 1 + z)/a 2 (x) dxdx' 
jff K w (x,x')K w (x,x")/o- 2 (x) dxdx'dx" 



(11) 



(12) 



where all integrals are in the range (0, z max ). The noise term has 
dropped out since its expectation value is zero. Adding Wq* t0 both 
sides of Eq. tl H we obtain 



fit 
w 



$ m (z')w tTue (z')dz' 



(13) 



The fitted constant, w *, is just a weighted average of the true equa- 
tion of state, with weighting function This weighting function 
is shown in Fig. 3 for the case of supernovae distributed evenly 
from z = to z = 2.0. It decreases steadily with redshift, as might 
be expected from the fact that all the supernovae are affected by 
the equation of state at z = 0, whereas the equation of state at 
z > 2.0 affects no supernovae. This result shows that the value of 
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Redshift z 

Figure 3. The weighting function for the constant w estimator for super- 
novae distributed uniformly up to z = 2.0. The solid line shows the weight- 
ing function computed with w = — 1 and the dashed one for w = 0. The 
small difference between these plots shows that the linear approximation 
holds to a good accuracy. 



wg obtained by fitting to supernovae data remains invariant if the 
w tTue (z') is changed in such a way that the weighted integral is 
unaffected. 

If the linear response approximation did not hold to a high 
accuracy then Eq. 1 1 31 would lead to different results depending 
on the value for w (id used to compute the weighting function $ w . 
Therefore a further test for the linear response approximation is 
provided by the variation of the weighting function evaluated for a 
different w hd . In Fig 3 we show the weighting function computed 
for two different values of w fld , namely w Rd = — 1.0 and w Rd = 0. 
The figure shows that the variation in the weighting function & w is 
small. 

To further quantify this difference we consider the change in 



w with the change in the weighting function 



Awq* = / A$ w {z')w ttue {z')dz 



(14) 



Squaring both sides and using the Schwartz inequality we obtain 
the bound: 

(Am?) 2 < /' ZmaX (A$,„) 2 dz [ Zn ™ X (w tIua ) 2 dz (15) 
Jo Jo 

Since we have assumed that \w\ < 1 in the relevant range, the 
second integral on the right hand side has a maximum value z max . 
(We have computed the weighting function for z max = 2.0.) With 
these numbers we can now compute (Auj fit ) 2 by considering the 



difference of evaluated at wo 



-1 and top — 0. We obtain 



\Aw b 



0.1. This shows that the expected error in the determina- 



tion of w through the linear response approximation is about 10% 
if we linearize the luminosity distance function around an incorrect 
initial guess. In practice, of course, one could do much better by 
first fitting the cosmological data with a constant w and then ex- 
panding around that value so that the approximation works better 
than when it is expanded about an arbitrary point. 



ever, quite cumbersome. We can simplify the notation a little by 
considering discretized expressions. Suppose the luminosity dis- 
tance is known at a large number of, uniformly distributed redshifts 
Zi where i — 1,N. We consider a model w(z) which is given at 
redshifts z' k where k = 1, M. 

After having chosen a w hd , which for simplicity is taken to be 
a constant, we can define the vectors d = {SDL(zi)} and w = 
{5w(z' i )} and matrix K, = {K(zi, z'j)Sz}, where /C is a N x M 
matrix and Sz denotes the redshift interval between the bins. With 
these definitions Eq {7) becomes 

d = JCw, (16) 
where both d and w are small quantities. Next we give the equiva- 
lent discrete version of the weighting function obtained in the pre- 
vious subsection. We have w Rt (z) = w Rd + Aw, where Aw is 
a constant; therefore, using the notation defined above we have 
w flt = uAw, where u = (1, 1, • • • , 1), and similarly w truc = 
w ttd + 5w(z), w truc = {Sw(zi)}. With these definitions Eq. {Tl} 
and 1 1 21 translate to 



U M true 

Aw = T w 
u 1 Mu 



A4 — IC W Id, 



If instead we fit the data with a linear model w 



(17) 



Awo+ Aw\z then, following the same procedure that led to Eq. I1 II 
we find 



Awo 



\ = 



[(z T Mz)u T M ~ (u T Mz)z T M] 

(u T 7Wu)(z T Xz) - (u T Mz) 2 
[{u T Mu)z T M - (u T Mz)u T M] 



(18) 



true / 1 ri\ 

w , (19) 



(u T Mu)(z T Mz) - (u T Mz) 2 

where we have also defined the vector z = {z^. This generalizes 
the concept of the weighting function to the linear case. It should 
be noted that the weighting function for Awo is different than in 
Eq. O. 

since we are also fitting for W\, However, as expected, 
these results agree for the special case when w truc = au. 

In the most general case we could approximate w(z) as a lin- 
ear combination of an arbitrary set of functions Fi (e.g. Gerke & 
Efstathiou, 2002) as 



fit/ \ fid 
W (Z) — W 



^CiFiiz) 



(20) 



In this case the coefficients c^s are related to the true equation of 
state as follows. Defining the N vectors Fi = {Fi(zk)} we obtain 
w flt = 53i=i c i^i- On performing a maximum likelihood analysis 
we find an expression for the coefficients a as follows 



where we have defined the two matrices 



= FlMF 
= FjjWw 



j 

true 



(21) 



(22) 
(23) 



4 EFFECTIVE w SEEN BY THE CMB 



3.2 Generalization 

We have derived above an expression for the expectation value of 
w m for the case of a constant w fit. The corresponding expres- 
sions for the case of more complicated fitting functions are, how- 



Current and near future supernova surveys will probe redshifts only 
up to around z = 2. We can also compute the effective equation 
of state probed by the Cosmic Microwave Background (CMB) at 
Zcmb ~ 1000. In the simplest case the CMB gives us the angular 
distance to the last scattering surface and therefore it does not give 
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Redshift z 

Figure 4. The weighting function which relates the effective equation of 
state ui c ff and the true equation of state w(z) as described in the text. The 
solid line shows the weighting function computed from the linear response 
approximation using u> fld = — 1 while the dotted one shows the weighting 
function inferred from £7q. The two functions agree to a good accuracy. 



1 2 5 10 20 

Redshift z 

Figure 5. The weighting function for the CMB evaluated using to fld = 
(solid line). The dotted line again shows the weighting function inferred 
from Qq. The two functions do not agree to a good accuracy, as discussed 
in the text. 



detailed information on the behaviour of w(z). Huey et al. (1999) 
show that a w c ff that keeps the CMB temperature and the matter 
power spectrum unaffected to an accuracy of 5-10%, is given by 
w ff = JilQ(a)w{a)da 

where the integral extends up to the CMB redshift. This is in the 
same spirit as the weighting function corresponding to fitting a 
constant constant equation of state described earlier. Following a 
similar procedure to that in the previous section we obtain 



W cft 



$ Cmb (2) 



$ culh (z)w(z)dz 



K'{l+z cmh ,l+z) 
jw, K'(l + z cmb ,l + z)dz 



(25) 



(26) 



where K' — K/(l + z) 2 due to the difference between Dl and 
D A . 

In Figure |4| (solid line) we plot this weighting function as 
calculated using the linear approximation by expanding around 
w — — 1. It has the same rough shape as the weighting function 
for supernovae but this time extends to higher redshift. We see that 
despite the high redshift of the CMB the weighting function falls 
off close to zero by redshift z ~ 5. This is due to the fact that the 
dark energy term in the Friedman equations is increasingly unim- 
portant as we go to high redshifts. For example it is well known 
that the cosmological constant is dynamically unimportant at CMB 
redshifts. 

Overlaid on Fig.|4](dashed line) is the weighting function from 
Eg. 1241 converting the integral from the scale factor a to redshift z. 
The good agreement shows that weighting with Qq works well for 
models close to the cosmological constant. In Fig [5] we show the 
same plot but this time expanding around w hd — 0. The two results 
do not agree to a good accuracy. This shows that the weighting with 
fig should not work so well for models which have an average 
w closer to zero. These figures also show that the linear response 
approximation works less well for the CMB since the weighting 
function changes by a large amount when we expand around w = 
— 1 and w = 0. This means that the weighting function would only 
apply for those models for which \w — w e ff] is small. 



5 NON-PARAMETRIC RECONSTRUCTION OF THE 
EQUATION OF STATE PARAMETER 

Although the fitting functions that we have considered so far pro- 
vide useful information about the true equation of state of the dark 
energy, the linear approximation discussed above in principle al- 
lows a non-parametric reconstruction of the equation of state. 

Consider a given SNe data set. To a first approximation we can 
fit it with a constant w, even though this may not be a good fit. We 
set the fiducial equation of state iu fld to this best fit constant value in 
all that follows. We next compute the difference between the SNe 
distances and the Dl computed for w fld to obtain the residuals 
SDl- We use the notation described in the first paragraph of § 13.21 
so these are given by d s , where the superscript s signifies that these 
are obtained from measured supernovae distances. These are re- 
lated to the w(z) in bins through Eg. 1161 where w = {w (zk) — wo} 
with wo given by the constant w fit. Since these are noisy estimates 
of distances we need to minimize the \ 2 function with respect to 
w to obtain the Maximum Likelihood estimator for w; that is, we 
minimize \ 2 = (d s — /Cw) T (d s — /Cw) and obtain 



{jc T icy 



(27) 



where we have assumed the noise on d" is constant with redshift, 
although the result could be extended for the general case. 

This equation also allows the trivial calculation of the Fisher 
matrix T = /C T /C corresponding to the uncertainties on the re- 
constructed equation of state. It was noted by Huterer & Starkman 
(2002) that the Fisher matrix was a surprisingly weak function of 
the model parameters which matches with the result obtained above 
and shows that it reflects on the validity of the linear response ap- 
proximation. 

Although Eq |27|g ives a formal solution of the problem, this es- 
timation is, in general, very noisy. Huterer et al. show that even for 
SNAP-like data there are only a few principal components which 
are well determined by this method. It is clear that this lack of res- 
olution is largely due to the fact that no constraints are imposed on 
the behaviour of w. We discuss some of the ways of rectifying this 
in a separate paper. 
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6 CONCLUSIONS 

We have shown that in the relevant range of parameters expected 
for the dark energy the luminosity distance is a linear functional 
to the equation of state to a surprisingly high level of accuracy. 
This approximation allows us to find the relationship between the 
usual polynomial models for w(z) and the true underlying equa- 
tion of state of the dark energy. Although the usual interpretation 
of these polynomial approximation is in terms of a series expansion 
of w(z), we show that the coefficients of the polynomial approx- 
imation are related to the true equation of state through an inte- 
gral relation. Only in the exact polynomial-like equation of state 
would these approximations measure the real w(z). We show that 
the fitting of cosmological data with such forms is still useful since 
the parameters of such models measure certain, well-defined, inte- 
grated properties of the underlying equation of state. Finally, the 
approximation allows a formal, non-parametric way to measure the 
equation of state. 
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